Explosive fragmentation of thin ceramic tube using pulsed power 
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This study experimentally examined the explosive fragmentation of thin ceramic tubes using 
pulsed power. A thin ceramic tube was threaded on a thin copper wire, and high voltage was applied 
to the wire using a pulsed power generator. This melted the wire and the resulting vapor put pressure 
on the ceramic tube, causing it to fragment. We examined the statistical properties of the fragment 
mass distribution. The cumulative fragment mass distribution obeyed the double exponential or 
power-law with exponential decay. Both distributions agreed well with the experimental data. We 
also found that the weighted mean fragment mass was scaled by the multiplicity. This result was 
similar to impact fragmentation, except for the crossover point. Finally, we obtained universal 
scaling for fragmentation, which is applicable to both impact and explosive fragmentation. 

PACS numbers: 46.50.+a, 62.20.Mk, 64.60.Ak 



Fragmentation is one of the most ubiquitous phenom- 
ena in nature 0. Well known examples include nuclear 
collision, the impact of brittle materials, and meteorites 
formed from asteroids. Fragmentation phenomena occur 
at all scales, from micro to macro, and most occur very 
rapidly. Consequently, it is difficult to observe fragmen- 
tation phenomena directly. The statistical properties of 
the fragment mass (or size) distribution have been inves- 
tigated frequently, and a scaling method is very useful for 
the analysis. This method is also used to study earth- 
quakes, turbulence, and many other non-equilibrium phe- 
nomena. Since the scaling laws are universal in many 
cases, we expect analogies across such phenomena. In 
particular, brittle fragmentation is a very simple exam- 
ple. Therefore, we believe that investigation of brittle 
fragmentation is indispensable for understanding non- 
equilibrium phenomena. 

Fragment mass distribution obeys the power-law form 
N(m) = J m n(m')dm' ~ m 1_r , that is, n(m) ~ m~ r , 
where m, n(m), and r are the fragment mass, fragment 
number of mass m, and characteristic exponent, respec- 
tively. Simple experiments examining brittle fragmen- 
tation have indicated the universality of fragmentation 
phenomena with this type of distribution 0, y, 0] ■ Sub- 
sequently, numerical simulations have been performed 
II B II 11 11 El El Recently, Astrom et al. proposed 
a generic scaling form, which consists of a combination 
of the power-law together with exponential decay and 
exponential functions [|J. This fits many classes of ex- 
periments and simulations, because it contains several 
typical function forms as limiting cases. 

Kun and Herrmann studied the transition between the 
damaged and fragmented states ||. This problem is the 
key to understanding the universality of fragmentation. 



They proposed that the transition involves the universal- 
ity of percolation. By contrast, Astrom et al. considered 
another critical behavior in the fragmentation transition 
0- Katsuragi et al. carried out experiments to check 
these predictions 0, 0. They also revealed the role 
of the log- normal distribution. A binomial multiplicative 
model was able to explain the log-normal distribution 
and multi-scaling nature of impact fragmentation. 

Through these simulations and experiments, the out- 
line of brittle fragmentation has been clarified systemat- 
ically; however, many unsolved problems remain. One 
that has received recent attention is the fragmentation 
of shells 0| . Wittel et al. examined impact and explo- 
sive fragmentation of eggshells in both simulations and 
experiments. They reported that explosive fragmenta- 
tion exhibits an abrupt transition, whereas impact frag- 
mentation shows a continuous transition. These results 
suggest that there is a new kind of universality that re- 
veals the explosive fragmentation of shells. In addition, 
a double exponential (not the power-law) fragment mass 
distribution has been observed in explosive fragmentation 
experiments |l 5| and simulations |lfj| . These results indi- 
cate the unresolved features of explosive fragmentation. 
Therefore, more investigations of explosive fragmentation 
are needed to consider fragmentation universality. 

In general, explosive fragmentation experiments are 
difficult to conduct because of the danger of the explo- 
sion. We performed relatively safe explosive fragmen- 
tation using a pulsed power generator. In this Letter, 
we report the results of our experiments. Moreover, we 
compare the results with previous impact fragmentation 
experiments. Finally, we discuss the universal scaling of 
impact and explosive fragmentation. 

We used a pulsed power generator consisting of four 
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FIG. 1: Schematic illustration of the experimental apparatus. 
A ceramic tube was threaded on a copper wire, and a pulsed 
high voltage was applied by the pulsed power generator. 
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capacitors with a capacitance of 0.1 /iF each, i.e., 0.4 /jF 
in total FigureQ]shows a schematic image of the ex- 
perimental apparatus. The capacitors were charged to a 
high negative voltage (around (— )20 kV), and triggered 
by an external signal to a gap switch. We fragmented 
thin 99.6% A1 2 3 ceramic tubes (SSA-S: Nikkato Cor- 
poration) that measured 1.2 mm in outside diameter, 0.2 
mm thick, and 200, 150, or 100 mm long, with a density 
of 3.9 g/cm 3 . They were sufficiently long that one might 
regard them as one-dimensional objects. In fact, a few 
(low discharge energy) cases produced one-dimensional- 
like fragments that retained the long tube geometry. 
However, we confirmed that most of the fragments had 
two-dimensional geometry, i.e., after fragmentation, very 
few fragments retained the original tube geometry, par- 
ticularly not those in the small mass range. Therefore, we 
considered this experiment a two-dimensional one. Ow- 
ing to the discharge voltage limit of the generator, we 
could not break larger samples than those used. 

A thin ceramic tube was threaded on a copper wire 
(0.05 mm diameter) and the wire was fixed to cop- 
per electrodes in a vessel. When a pulsed electric dis- 
charge passed through the wire, the copper vaporized 
and the resulting pressure inside the ceramic tube in- 
creased dramatically, causing fragmentation of the ce- 
ramic tube. During the experiments, the current and 
voltage between the electrodes were measured using a 
current probe (Model 110A: Pearson Electronics) and 
a high-voltage probe (EP-50K: Pulse Electric Engineer- 
ing), respectively. After each fragmentation, all of the 
fragments were collected and their masses were measured 
using an electronic balance (ER-60A: A&D). We frag- 
mented 47 ceramic tubes (21: 200 mm-long, 18: 150 
mm-long, and 8: 100 mm-long). The input parame- 
ter of this experiment was the discharge voltage, which 
ranged from 18 ~ 23 kV. When we analyzed the frag- 
ment statistics, we used the minimum limit of fragment 
mass TOmin = 0.001 g. This roughly corresponds to 
the reliable limit of our measurements. Moreover, this 
value is sufficiently larger than the smallest limit mass 
of two-dimensionality 0.2 3 (mm 3 ) x3.9 x 10~ 3 (g/mm 3 ) 



FIG. 2: Typical profiles of the current I, voltage V, power P, 
and discharge energy E as functions of time. A sharp power 
pulse can be confirmed on an approximately /is time scale. 

= 3.1 x 10~ 5 g. 

First, we show typical profiles of the current, voltage, 
power, and released energy during a discharge in Fig. 
Other experiments showed similar profiles qualitatively. 
A very sharp power pulse was confirmed. The electrical 
discharge happened on a fis time scale, which resulted 
from the electrical properties of the circuit. This resulted 
in efficient melting of the copper wire. 

Next, we show merged cumulative fragment mass dis- 
tributions in Fig. |3 In Fig. the N(m) curves are de- 
fined as N(m) — J2l=i ^i{ m )l r i where r and Ni(m) are 
the number of sum-ups and the zth cumulative distribu- 
tion, respectively, i.e., the N(m) curves are averaged. We 
averaged the results under the same conditions (discharge 
voltage and sample size). The mass to was normalized 
to the average total mass mtot ■ For the cases with a high 
voltage discharge, the N(m) curves show a power-law-like 
portion. This regime corresponds to the fully fragmented 
state. By contrast, with low voltage discharge, the N(m) 
curves appear almost flat in the small mass region. They 
clearly differ from the fully fragmented regime distribu- 
tions, and are similar to an integrated log-normal dis- 
tribution. These features resemble those of impact frag- 
mentation 0, 0| . We also show the relation between 
the estimated discharge energy per unit mass (Ed/rritot) 
and the average total fragment number per unit mass 
(N/mtot) in the inset of Fig. El Roughly speaking, the 
estimated discharge energy shows the extent of the in- 
ner pressure produced, and the average total fragment 
number indicates the degree of fragmentation. 

We tried to fit the data in the fully fragmented regime 
using two traditional forms: the double exponential 
form N(m) = Ai exp(— m/mi) -I- A2 exp(— 111/7112), and 
the power-law with exponential decay form N(m) — 
Aom^ 1 " -1 ) exp(— m/mo). Where r, too, mi, to-2, Aq, 
A±, and A2 are fitting parameters. The results of fitting 
for the 150-mm-long and 21-kV-discharge case are shown 
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FIG. 3: Typical cumulative fragment mass distributions. The 
distribution curves were obtained from a mean of 4, 5, 3, and 6 
fragmentation events for 20 kV (150 mm), 21 kV (150 mm), 20 
kV (200 mm), and 23 kV (200 mm), respectively. The inset 
shows the relation between the estimated discharge energy 
and average total fragment number. 

in Fig. 0] Both function forms agree with the experimen- 
tal results very well. Other N(m) curves for the fully 
fragmented regime were also fitted using both forms. We 
cannot conclude which form is more appropriate to ex- 
plain the experiment. 

The generic scaling form proposed by Astrom et al. 
(Eq. (3) in Ref. Q) is close to the double exponential 
form when the penetration depth in the scaling form 
is small. The penetration depth indicates the strength 
of crack progression whose dynamics result in a power- 
law fragment mass distribution In our experiment, a 
small penetration depth is reasonable because the pres- 
sure loads were applied uniformly to the two-dimensional 
(2D) space of the samples. Fragmentation occurred si- 
multaneously throughout 2D space due to this loading 
condition symmetry. Then, crack nucleation determines 
the fragment mass distribution rather than the deep 
crack dynamics. Consequently, the fragment mass distri- 
bution approximates the double exponential form due to 
the Poisson random nucleation of a crack and the dynam- 
ics of neighboring cracks that lack a sufficient power-law 
(i.e., a small penetration depth). In fact, the horizontal 
sandwich experimental results in Ref. [12j can be fitted 
using Astrom et al.'s scaling form when that is close to 
the double exponential form [l^. The horizontal sand- 
wich experiment has similar loading condition symmetry 
to the present explosive experiment. 

Conversely, an empirical form of the power-law with 
exponential decay also results in good agreement. This 
form has fewer fitting parameters than the double expo- 
nential form and is consistent with other experiments, 
analyses, and simulations performed to date. A weak- 
ness of this fitting is the narrow region of the power-law 
regime in our data. Since the linear part of the log-log 
plot covers less than one order of magnitude, it might not 




FIG. 4: Fitting results using the double exponential and 
power-law with exponential decay forms. Circles indicate the 
result for 150-mm-long tubes with a 21 kV discharge. Both 
fittings make good agreement. 



be sufficient to conclude that the power-law form is valid. 
Note that the double exponential function can produce a 
similar portion, as seen in Fig. 0] We found r — 1 values 
in the range 0.66 — 0.76 for the fully fragmented regime. 
These values differ from the results of Wittel et al. [lij . 
There are two free boundaries in the geometry of a tube, 
whereas an eggshell has only one. To check the boundary 
effect, we also experimented with two-tone colored tubes. 
However, we found no significant difference between the 
regions near the ends and that in the middle of the tubes. 
While Wittel et al.'s explosion arose from a point source 
and expanded spherically, we used a linear source and 
the explosion expanded cylindrically. The curvature of 
the samples might also have affected the results. These 
reasons might explain the differences in the results, al- 
though we think that the problem remains unresolved. 

We also examined the behavior of the fcth order mo- 
ment Mk = Y2m Tn k n(m) in order to elucidate the criti- 
cal scaling in the explosive fragmentation. Note that this 
summation involves the lar gest fragment mass. Similar 
to the previous analyses 0,^3 > we show the log-log plot 
of M* = 7\f 2 /(M 1 m min ) vs. fi = M m min /Afi in Fig. 
Figure also shows previous results for impact fragmen- 
tation [121 Il3| , where M* and fj, are the wei ghte d mean 
fragment mass and multiplicity, respectively Jl^L Il3l Hif . 
When we assume the scaling ~ yT a , we can confirm 
the crossover of the scaling exponent a at a certain fi c 
[lflj . As Fig. El shows, a crossover from a small a value 
to a large one can also occur in explosive fragmentation. 
Nevertheless, the crossover point fi c does not seem to be 
universal for impact and explosive fragmentation. Al- 
though we have conjectured on the universality of /i c in 
Ref. [l3j, it shows non- universal behavior. However, this 
is not surprising. The slopes (exponent) a appear to be 
universal in the two. We think that this is reasonable 
in terms of critical universality. In general, the critical 
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FIG. 5: Log-log plot of the weighted mean fragment mass M* 
vs. the multiplicity fi. The crossover point \i c is not universal 
for impact and explosive fragmentation. However, the slopes 
(exponent) appear to be universal. 
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FIG. 6: Universal scaling plots of the weighted mean fragment 
mass. We use the values a = 0.4 and /3 = 0.7. We show (a) 
M*/L a vs. fi/L , and (b) M*/L a+0S as a function of p. 



point is not universal despite the universality of critical 
exponents. Moreover, this difference in fi c might arise 
from differences in the materials, i.e., glass or ceramic. 

We would like to determine a more universal law for 
fragmentation criticality. In order to improve the scal- 
ing, we use the normalized largest fragment mass L = 
m-h/mmim where mj, is the largest fragment mass. In the 
finite size system, L plays an important role for scaling 
analyses. Finally, we assume the universal scaling of the 
weighted mean fragment mass, as shown in Figs. [SI which 
include both impact and explosive results. The scaling 



form should be written as M*/L a — /(/i/L' 3 ). In Fig. 
Et a ), we use the values a = 0.4 and (3 = 0.7. From Fig. 
Eta), we deduce the scaling function f(x) ~ x~ s . We ob- 
tain the exponent 6 = 0.424 ± 0.006 from fitting the data 
in Fig.Et a )- Furthermore, this implies M*/L a+)3S ~ fi~ s . 
We show the log- log plot of Al r /L a+ ^ s as a function of [x 
in Fig.Etb). The 6 value from Fig.Etb) is S = 0.42±0.01. 
In either case, the scaling appears convincing globally. 
We believe that must be a crucial scaling for critical 
fragmentation. Specifically, the weighted mean fragment 
mass M* can be scaled as i°- 69 ^-°- 42 independent of the 
fragment mass distribution type: log-normal, power-law, 
or double exponential. The scaling is very robust for 
small to large imparted energies. 

If we assume the binomial multiplicative model with 
parameter a = 2/3 12] for the small fi regime, we can 
calculate <5 bin as 



a log a — \og[a 2 + (1 — a) 2 ] 



(1) 



log 2-/3 log a 

Substituting a — 0.4, (3 — 0.7, and a = 2/3 into 
Eq. we obtain <5 bin = 0.44. This value is close to 
5 = 0.42. By contrast, if we assume the power-law 
form (with 1 < r < 2) for the large fi regime, Mk can 
be expressed as ~ f n L m k mT T dm. From this no- 
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tation, we can approximate M*/L a 

[L(3-r)/£(2-r)]/[ L a] _ ([£(r-2)]/[ i /9] 



(^/L^rw i 

ipow. therefore, 



4,™ = 



1 - a 
2 + /?- 
0.4, 



(2) 



When we substitute a = 0.4, (3 = 0.7, and r = 1.7 
into Eq. (J2J, we are led to <5 pow = 0.6. This value is 
slightly different from S = 0.42. This suggests the limit 
of the power-law approximation, or the importance of the 
exponential form. Determining the origin of the values a, 
/3, and S remains a future problem. A theoretical study 
of this scaling would clarify the fragmentation criticality 
well. 

In conclusion, we carried out a two-dimensional explo- 
sive fragmentation experiment using pulsed power. The 
N(m) curves for the fully fragmented regime could be fit- 
ted using the double exponential form or the power-law 
with exponential decay. Each form has merits and de- 
merits. The weighted mean fragment mass M* was scaled 
using the multiplicity /i and its scaling showed crossover, 
just as in impact fragmentation. While the scaling ex- 
ponent a was the same as for impact fragmentation, the 
crossover point /i c differed. Using the normalized largest 
fragment mass L, we obtained universal scaling between 
the impact (both vertically and horizontally) and explo- 
sive fragmentation data. This scaling also appears inde- 
pendent of the material broken. 
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